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Abstract. A pebbling move on a weighted graph removes some pebbles at a vertex and 
adds one pebble at an adjacent vertex. The number of pebbles removed is the weight 
of the edge connecting the vertices. A vertex is reachable from a pebble distribution 
if it is possible to move a pebble to that vertex using pebbling moves. The pebbling 
' number of a weighted graph is the smallest number m needed to guarantee that any 

^ . vertex is reachable from any pebble distribution of m pebbles. Regular pebbling problems 

^ I on unweighted graphs are special cases when the weight on every edge is 2. A regular 

pebbling problem often simplifies to a pebbling problem on a simpler weighted graph. 
We present an algorithm to find the pebbling number of weighted graphs. We use this 
' algorithm together with graph simplifications to find the regular pebbling number of all 

connected graphs with at most nine vertices. 
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! 1. Introduction 

> . 

j3 ' Graph pebbling has its origin in number theory. It is a model for the transportation of 

resources. Starting with a pebble distribution on the vertices of a simple connected graph, a 
pebbling move removes two pebbles from a vertex and adds one pebble at an adjacent vertex. 
We can think of the pebbles as fuel containers. Then the loss of the pebble during a move 
is the cost of transportation. A vertex is called reachable if a pebble can be moved to that 
vertex using pebbling moves. The pebbling number of a graph is the minimum number of 
VjD . pebbles that guarantees that every vertex is reachable. There are many different variations 

of pebbling. For a comprehensive list of references for the extensive literature see the survey 
tJ- ■ papers [H[5]. 

^ , Our goal is to find an algorithm that finds the pebbling number in a realistic amount 

of computing time. The main idea of the algorithm is that if we know all the sufficient 
distributions from which a given goal vertex is reachable then we can find the insufficient 
distributions from which the goal vertex is not reachable. An insufficient distribution must 
be smaller than a sufficient distribution. The pebbling number can be found by finding an 
. insufficient distribution with the most pebbles. The problem is that there are too many 

sufHcient distributions. Luckily it suffices to find the barely sufficient distributions from 
which the goal vertex is no longer reachable after the removal of any pebble. 

Our algorithm works even if the cost of moving a pebble from one vertex to another 
varies between different vertices. To take advantage of this, we generalize the notion of 
graph pebbling on weighted graphs and we develop the basic theory of weighted graph 
pebbling. 

The generalization is worth the effort since pebbling on many graphs can be simplified 
if we replace the graph by a weighted graph with fewer edges. For example a tree can be 
replaced by a weighted graph containing a single edge. Cut vertices, leaves and ears offer 
the most fruitful simplifications. 
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We use these simplifications and our algorithm to calculate the pebbling number of all 
connected graphs with at most nine vertices. We present the spectrum of pebbling numbers 
in terms of the number of vertices in the graph. 

2. Preliminaries 

Let G be a simple connected graph. We use the notation V{G) for the vertex set and 
E{G) for the edge set. We use the standard notation vu = uv for the edge {v,u} G E{G). 
A path of G is a subgraph isomorphic to the path graph P„ with n > 1 vertices. A weighted 
graph Guj is a graph G with a weight function uj : E{G) N. 

A pebble function on G is a function p : V{G) Z where p{v) is the number of pebbles 
placed at v. A pebble distribution is a nonnegative pebble function. The size of a pebble 
distribution p is the total number of pebbles \\p\\ = ^y^v(G)Pi''^)- '^^^ support of the 
pebble distribution p is the set supp(p) = {v G V{G) \ p{v) > 0}. We are going to use the 
notation p{vi, . . . , Vn, *) = {a-i, . . . , an, q{*)) to indicate that p{vi) = Oj for z G {1, . . . , n} 
and p{w) = q{w) for all w G V{G) \ {vi, . . . ,Vn}. 

If vu € E{G) then the pebbling move (v^u) on the weighted graph G^ removes lo{vu) 
pebbles at vertex v and adds one pebble at vertex u, more precisely, it replaces the pebble 
function p with the pebble function 

P{v^u){v,U,*) = {P{v) -Uj{vu),p{u) + l,p{*)). 

Note that the resulting pebble function P(v^u) might not be a pebble distribution even if p 
is. 

The inverse of the pebbling move (v^u) is denoted by {v-^u)~^ . The inverse removes 
a pebble from u and adds two pebbles at v, that is, it creates the new distribution 
(f , tt, *) = {p{v) + 2,p{u) — l,p{*)). Note that {v^u)~^ is not a pebbling move. 

A pebbling sequence is a finite sequence s = (si, . . . , s^) of pebbling moves. The pebble 
function gotten from the pebble function p after applying the moves in s is denoted by pg. 
The concatenation of the pebbling sequences r = (ri, . . . , rfc) and s = (si, . . . , s^) is denoted 
by rs = (ri, . . .,rk,si, . . .,si). 

A pebbling sequence (si, . . . , Sn) is executable from the pebble distribution p P(si,...,si) is 
nonnegative for alH G {1, . . . , n}. A vertex x of G is t-reachable from the pebble distribution 
p if there is an executable pebbling sequence s such that Ps{x) > 1. We say x is reachable 
if it is 1-reachable. 

We write 7rt{Gi^,x) for the minimum number m such that x is t-reachable from every 
pebble distribution of size m. We use the notation 7r(Ga,,x) for tti{Gu!,x). The t-pebbling 
number 7rt{Gt^) is max{7rt(Ga;, x) \ x E V{G)}. The pebbling number Tr{Gij) is the 1-pebbling 
number 7ri(Gt^,). 

If uj{e) = 2 for all e G E{G) then 7r{G^) = it{G) is the usual unweighted pebbling number. 
So we allow the weight function uj to be defined only on a subset of V{G) and use the default 
weight of 2 for edges where lo is undefined. 

Changing the order of moves in an executable pebbling sequence s may result in a sequence 
r that is no longer executable. On the other hand the ordering of the moves has no effect 
on the resulting pebble function, that is, p^ = Pr- This motivates the following definition. 

Given a multiset S of pebbling moves on the weighted graph (G^j), the transition digraph 
T{G,S) is a directed multigraph whose vertex set is V{G), and each move (v^u) in S 
is represented by a directed edge {v, u). The transition digraph of a pebbling sequence 
s = (si, . . . , Sn) is T(G, s) = T{G, S), where S = {si, . . . , s„} is the multiset of moves in 
s. Let d~^^ denote the in-degree and d^^^ ^-^ the out-degree in T(G, S). We simply write 

and d^ if the transition digraph is clear from context. It is easy to see that the pebble 
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function gotten from p after applying the moves in a multiset S of pebbling moves in any 
order satisfies 

ps{v) =p{v)+d-^^^^^{v)-Y,Mv^) I {v,u)eE{T{G,S))} 
for all V & G. For unweighted graphs the formula simplifies to 

Ps{v) = p{v) + d'^^^g^iv) - 2d+(g__5-)(t;). 

3. Cycles in the transition digraph 

In this section we present a version of the No-Cycle Lemma [U El [9]. If the pebbling 
sequence s is executable from a pebble distribution p then we clearly must have ps > 0. We 
say that a multiset S of pebbling moves is balanced with a pebble distribution p at vertex v 
if ps{v) > 0. The multiset 5 is balanced with p if 5 is balanced with p at all v £ V{G), that 
is, PS > 0. We say that a pebbling sequence s is balanced with p if the multiset of moves 
in s is balanced with p. The balance condition is necessary but not sufficient for a pebbling 
sequence to be executable. A multiset of pebbling moves or a pebbling sequence is called 
acyclic if the corresponding transition digraph has no directed cycles. 

Proposition 3.1. If S is a multiset of pebbling moves on G^^ then there is an acyclic multiset 
i? C S such that pR > ps for all pebble function p on G. 

Proof. Let p be a pebble function on G. Suppose that T(G, S) has a directed cycle C . Let Q 
be the multiset of pebbling moves corresponding to the arrows of C and R = S\Q. Let 
be the first vertex from v along G . Then pr{v) = psiv) — 1 + u){vuv) > ps{v) for v € V{G) 
and pr{v) = ps{v) for v G \ V{G). 

We can repeat this process until we eliminate all the cycles. We finish in finitely many 
steps since every step decreases the number of pebbling moves. □ 

Definition 3.2. Let S* be a multiset of pebbling moves on G. An element (v^u) S 5" is 
called an initial move of 5 if d~{v) = 0. A pebbling sequence s is called regular if Si is an 
initial move of 5 \ {si, . . . , Si-i} for all i. 

It is clear that if the multiset S is balanced with a pebble distribution p and s is an initial 
move of S then s is executable from p. 

Proposition 3.3. IfS is an acyclic multiset then there is a regular sequence s of the elements 
of S. If S is also balanced with the pebble function p then s is executable from p. 

Proof. If S is acyclic then we must have an initial move t of S. Then S \ {t} is still acyclic. 
So we can recursively find the elements of s recursively by picking an initial move t of S and 
then replacing S with 5* \ {t} at each step. 

Now assume that S is balanced with p. Then Si = S \ {si, . . . , Si-i} is balanced with 
for all i since {p(si,...,Si--i))Si = > 0. Hence the initial move Sj of Si is executable 
from that is, P{si,...,s,) > for all i. □ 

The following result is our main tool. 

Theorem 3.4. Let p be a pebble distribution on G^^ and x e y{G). The following are 
equivalent. 

(1) Vertex x is reachable from p. 

(2) There is a multiset S of pebbling moves with P5 > and ps{x) > 1. 

(3) There is an acyclic multiset R of pebbling moves with pn > and pr{x) > 1. 

(4) Vertex x is reachable from p through a regular pebbling sequence. 
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H „ K 

Figure 4.1. Simplification using the cut vertex v. If 7r((K^, v) = at + 6 then 

'n{Gui,x) = tt{Gcj,x) + b. 

Proof. If X is reachable from p then there is a sequence s of pebbling moves such that s is 
executable from p and ^^(a;) > 1. If 5 is the multiset of the moves of s then ps > and 
Ps{x) > 1 and so (1) implies (2). 

By Proposition [3TT], (2) implies (3) and by Proposition 13.31 (3) implies (4). It is clear that 
(4) imphes (1). □ 

It is convenient to write the condition ps > and ps{x) > 1 compactly as ps > l{x} 
using the indicator function of the singleton set {x}. 

4. Cut vertices 

The pebbling number of a graph with a cut vertex often can be calculated using a simpler 
graph. This simplification introduces new weights. The following theorem is the main reason 
we study weighted graphs. 

Proposition 4.1. Let H and K he connected graphs such that V{H)r\V{K) = {v} and v is a 

cut vertex ofG = HUK. Letuj he a weight function onE{G). Assume that TTt{K^,v) = at+b 
for all t. Define a graph G hy V{G) = V{H)U{u} and E{G) = E{H) U {vu}. Define a 
weight function on E{G) hy 

. , , I a if e = vu 
lLi;(ej else 

If the goal vertex x is in V{H) then tt{Guj,x) = t:{Gcj^x) + b. 

To simplify notation, we used instead of the more precise -ftr^|_E;(ii-) even though uj is 
defined on values outside of E{K). 

Proof. The graphs are visualized in Figure [411 First we show that 7r((S'^,a;) > -KiGcb^x) + b. 
Let p be a pebble distribution on G with \\p\\ = 7r(G^,x) — b. We create a new distribution 
g on G consisting of red and green pebbles. The red pebbles are placed on H exactly the 
same way as the pebbles in p are placed on H . The number of green pebbles is p{u) + b. 
The green pebbles are placed on K so that the number of pebbles that can be moved to v 
using only green pebbles is minimum. This minimum number is clearly lp{u)/a\. Note that 
q can have both red and green pebbles on v. Then ||g|| = 7r(G(^,x), so there is an acyclic 
multiset S of pebbling moves on G^j such that qs > l{x} ■ Let Sh and Sk contain the moves 
of S inside H and K respectively so that S = Sh^Sk- 

We are going to see that moving red pebbles from H to K is not beneficial and so these 
moves can be eliminated. Since S is acyclic, any maximal walk in T{K,Sk) starting at v 
is actually a path. Let us remove the pebbling moves corresponding to such maximal walks 
from Sk until we eliminate all walks from T{K,Sk) starting at v. The choice of these 
maximal walks is not unique and they can overlap, so we need to eliminate them one by 
one. The resulting multiset S'^ C Sk is balanced with q and qsxi'^) ^ 95' (^)- 

Executing S*^ from q cannot move more than \_p{u)/a\ pebbles to v since 5^ does not 
have any effect on the red pebbles. Let R be the multiset containing the elements of Sh 
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together with LIIp(^) ll/f^^J copies of the move Then it is clear that pr{u) > and pR 

is not smaller than qs on H, hence pR > l{x}- 

Now we show that 7r{Guj,x) < ■K{Gcj,x) + b. Let p be a pebble distribution on G with 
IIpII = TT{Gcj,x) + b. Let c = be the number of pebbles in this distribution on 

V{K) \ {v}. We create a new distribution g on G such that q and p are the same on V{H) 
and q{u) = max{0, c — 6}. Then \\q\\ > n^Gcjjx) so there is an acyclic multiset S of pebbling 
moves on G;;j such that qs > l{x}- There is a multiset Ri of pebbling moves on K using 
only the pebbles on V{K) \ {v} such that PR^iy) = p{v) + [max{0,c — b}/a\. Let R be 
the multiset containing the elements of Ri together with the elements of S different from 
{u^v). Then pR is not smaller than qs on H and pR is nonnegative on V{K) \ {u}, hence 

PR>l{a;}- □ 

Note 4.2. The previous proposition is applicable in many situations since the function t i— > 
7rf(G, x) is often linear; for example for trees, complete graphs and hypercubes. In particular 

it is hnear for cycles [6] where TTt{C2n) = t2" and ■nt{C2n+i) = I + {t - 1)2" + 2 ^ . 

The simplest nonlinear example is the wheel graph with 5 vertices. If x is a degree 3 
vertex then 

"*^^^'") = \4t if^>2- 
If G is the complete graph with 7 vertices with one missing edge xy then 

/2t + 5 ifiG{l,2} 
^ ' ' \At ifi>3 

5. Simplifications using leaves 

Proposition 5.1. Let G^ he a weighted star graph with center x and spikes xvi, . . . ,xVn- 
Suppose that a = ui{xvi) is the maximum value ofuj. Then 7rt{Guj,x) = ta+^^^2i^i^''^i)~^) ■ 

Proof. The maximum number of pebbles we can place on Vi so that at most ti pebbles can 
be moved from Vi to x is (ij + l)uj{xvi) — 1. So 

n 

MGu,,x) = max{^((ti + l)a;(xt;i) - 1) | ti + ■ ■ ■ + t„ < t} + 1 

1=1 



{t-l + l)uj{xvi) - 1 + ^((0 + l)io{xv,) - 1) + 1 



i=2 

n 

= ta + - 1) 

1=2 

since the maximum is taken when ti = t — 1 and t2 = • • • = tn = 0. □ 

The reader can easily verify the following result. 

Proposition 5.2. Let x, vi and V2 be the consecutive vertices of the graph G = P3 with 
weight function u. Then ■Kt{Guj,x) = tuj{xvi)uj{viV2) . 

Propositions 14. 11 15. II and [5?2] allow us to calculate the pebbling number of every weighted 
tree since we can simplify the tree to a single edge. The process is shown in the next example. 

Example 5.3. Figure ISTTl shows the stages of the simplification of a tree. First we let K be 
the subgraph of G generated by {f3,f4,f5}. Then irtiK^v^) = At by Proposition 15.21 so we 
replace K by the weighted edge ^3^1 to get G^V ■ Next we let K be the subgraph of G^V 
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ivo) IVo) 

VI—V3 V4 V5 V— V3 — Ui W V3 — U2 [Vo)—U3 

V2 ^ V2 ^ 



Figure 5.1. Simplification of a tree: 7r(G, wq) = vr(GL\\i'o) = 7r(G'i^2\^o) + 

.(3) 



2 = 7r(GS,r;o) + 2 = 10. 





Figure 6.1. Substitution for a closed ear. The edge weight a is 2^2 J where 
k is the number of vertices of the closed ear. 



generated by {vi,V2,ui} . Then TTtiK^v^) = 4t + 2 by Proposition 15.11 so we replace K by 

-.(3) 

6. Simplification using ears 



(o\ 

the weighted edge f 3ii2 to get Gu)2 ■ Finally we use Proposition 15.21 again to get Gl_ 



In this section we use the existence of special paths in our graph to simplify the calculation 
of the pebbling number. A thread of a graph is a path containing vertices of degree 2. 

Definition 6.1. Let x be a goal vertex in G. Let wi, . . . , f„ be the consecutive vertices of a 
maximal thread T not containing x. There are unique vertices vq and f^+i outside of T that 
are adjacent to vi and Vn respectively. The subgraph E induced by vq, • • • , Vn+i is called an 
ear. The vertices of T are called the inner vertices of E. If vq = Vn+i then E is called a 
closed ear. If the vertices of T are cut vertices then E is called a cut ear. If E is neither a 
closed ear nor a cut ear then it is called an open ear. 

Note that an ear has at least two edges. Also note that the goal vertex can be an end 
vertex of an ear. 

6.1. Closed ears. If a closed ear has default weights then it can be replaced by a weighted 
edge using Corollary 14.11 and Note 14.21 The simplification is shown in Figure 16.11 If the 
closed ear has 2n vertices then -KiG^^x) = TT{Gcj,x). If the closed ear has 2n + 1 vertices 

then 7r{Guj, x) = TT{Gcj,x) + 1 — 2" + 2 . The edge weight is a = 2" in both cases. 

6.2. Cut ears. Cut ears can be replaced by weighted edges as well. First we need the 
following result. 

Lemma 6.2. Let p be a maximum size pebble distribution from which the goal vertex x is 
not reachable. Then p has no pebbles on the inner vertices of a cut ear. 

Proof. Suppose u is an inner vertex of the cut ear E and p{u) > 0. Let H and K be the 
connected components of G \ {u} such that x £ H. There is a unique vertex v € K that is 
adjacent to u. The size of q = is larger than the size of p. We show that x is not 

reachable from q which is a contradiction. 

Suppose there is an acyclic multiset 5 of pebbling moves with qs > f{x}- If (v^u) G S 
then with R = S\{{v^u)} we have Pr = qs ^ which is not possible. So we can assume 
that ^ S. Let R contain those moves of S that do not involve any vertex in K. Then 
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Figure 6.2. Substitution for a cut ear. The edge weight is u;{vu) = a = 2"~^ 
where n is the number vertices of the path connecting w to n in G. 



V 
X 



Figure 6.3. A graph with an open ear. 

PRiu) > qs{u) + 1, Pr{w) = qs{w) for w £ V{H) and Pr{w) = p{w) for w G V{K). So 
PR > l{a;} which is again impossible. □ 

Proposition 6.3. Let E he a cut ear of G^^ with end vertices v and u. Let G he the graph 
created from G hy removing the inner vertices of E and adding the edge vu. Define Cj on 
E{G) hy 

a if e = vu 
oj{e) else 

where a is the product of the weights of the edges of E. If the goal vertex x is not an inner 
vertex of E then 7r(Gt^, a;) = 7r{Gcd,x). 



w(e) 



Proof. Without loss of generality we can assume that x is closer to v than to u as shown in 
Figure [6?2l Let u = vi,V2, ■ ■ ■ ,Vk = v he the consecutive vertices of E. 

First we show that Tr{Gi^,x) < tt{Gcj,x). For a contradiction, assume that 7r(G^,j;) > 
7r(G^, x). Let p be a maximum size pebble distribution on G from which x is not reachable. 
By Lemma [621 P has no pebbles on the inner vertices of E so the restriction q = is 

a pebble distribution on G with ||g|| = \\p\\ = 7r(G^,a;) — 1 > Tr{Gcj,x). Hence there is a 
multiset S of pebbling moves on Gcj such that qs > Let R be the multiset of pebbling 
moves containing the moves in S with each move of the form (n->f ) replaced by the moves 
{vi^V2), . . . , (vk-i^Vk)- Then qji > l^^y which is a contradiction. 

Now we show that 7r(Gi^,x) > TT{Gcj,x). Let p be a pebble distribution on G with size 
ir{Goj,x). Let q be the extension of p to V{G) such that q is zero on the inner vertices of 
E. Then ||g|| = = it{Guj,x) and so there is an acyclic multiset S of pebbling moves such 
that qs > We create a multiset R of pebbling moves on G as follows. We start with 

5. We then search for a directed path in T(G, S) connecting u to v and we remove all the 
moves corresponding to the arrows of this directed path. We do this until there are no more 
such directed paths. Then we add as many copies of as the number of directed paths 

removed. Finally we remove all moves involving inner vertices of E. It is easy to see that 
PR>'i-{x}- □ 



6.3. Open ears. Figure 16.31 depicts an open ear. An open ear cannot be replaced by a 
single edge but we can still take advantage of it using squishing as explained in Section [H 

6.4. Examples. Figures [6^ and [6?5] show two examples of simplified graphs using ears. The 
graph G is the same in both examples but the ears are different because the goal vertices 
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(fOh -Vi- -V2 
f 5 - - 174 



^^3 ::: 



, Vg ::::::::: Vj 



- Vg ::::::::: ^8 



G 



®< 4 

V3 U 

G,7, 



Figure 6.4. Simplification of G with goal vertex vq such that tt{G,vo) = 
TT{Gcj,vo) + 1 = 36. The two cut ears denoted by dashed edges are replaced 
by the weighted edges vqv^ and ^5^3 in G. The closed ear denoted by double 
dotted edges is replaced by the weighted edge v^u. 



Vo - -Vi- -V2, 
V5- -V4 



V3 



ive) V7 

■ Vg V8 



V3 



V5 4 



iV6) V7 

■ Vg V8 



G 



Ga 



Figure 6.5. Simplification of G with goal vertex fg such that 7r(G, fg) = 
'k{Gcj-,v%) = 22. The two cut ears denoted by dashed edges are replaced by 
the weighted edges vqV3 and v^V3 in G. The open ear denoted by dotted 
edges remains in G. 



are different. The pebbling number of the graph is 7r(G) = 36. The path connecting ^3 to 
could have been simplified using leaves in both examples. The path connecting vq to ^3 

could have been simplified using leaves as well, but only in the second example. In both of 

these cases a further simplification is possible using leaves and Proposition 15.11 

Note that in the second example the end vertices ^3 and of the open ear are adjacent. 

This possibility is important to keep in mind during the development of an algorithm to 

find open ears. 

7. Barely sufficient pebble distributions 

Let T>{G) be the set of pebble distributions on the graph G. For p,q £ T^{G) we write 
p < q if p{v) < q{v) for all v G G. This gives a partial order on T>{G). We write p < q if 
p < q but p ^ q. It is clear that if a goal vertex is reachable from p and p < q then the goal 
vertex is also reachable from q. 

Definition 7.1. Let x be a goal vertex of Goj- A pebble distribution p is sufficient for x if a; is 
reachable from p. The set of sufficient distributions for x is denoted by S{Guj,x). A pebble 
distribution p G S{Guj,x) is barely sufficient for x if x is not reachable from any pebble 
distribution q satisfying q < p- The set of barely sufficient distributions for x is denoted by 
x). The set of insufficient distributions for x is T{G^,x) = ^{G) \ S{Gu}, x). We are 
going to use the notation 5(x), B{x) and I{x) if G and to is clear from the context. 

We can partition B{x) into the disjoint union l3o{x)0 ■ ■ ■ iJBk{x) where Bi{x) contains 
those distributions in B{x) from which x is reachable in i pebbling moves but x is not reach- 
able in fewer than i moves. Note that the only element of Bo{x) is the pebble distribution 
that contains a single pebble on x. 

Example 7.2. Figure ET] shows an example of B{G,x). 

The following result is the main reason for our interest in barely sufficient distributions. 

Proposition 7.3. We have p G T(G^,x) if and only if q < p does not hold for any q G 
BiG^,x). 
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G l3o{G,x) Bi{G,x) B2{G,x) Bs{G,x) 

Figure 7.1. The barely sufficient pebble distributions for vertex x. The 
vertices denoted by bullets have no pebbles. 

Proof. If X is reachable from a pebble distribution p then we can remove pebbles from p one 
by one if needed until we get a g G B{G^,x) that satisfies q < p. The other direction of the 
result is obviously true. □ 

The following example shows how Proposition 17.31 can be used to find the insufficient 
distributions. 

Example 7.4. In Example 17.21 the maximal elements of I{G,x) are p{x,vi,V2,V3) = 
(0, 1, 1, 1), q{x, fi, f2, fs) = (0, 0, 3, 1) and r{x, vi,V2,V3) = (0, 0, 1, 3). The maximum size is 
|g| = 4 = |r| and so tt{G,x) = 5. 

Our purpose now is to construct algorithms for finding B{G^^,x) and 7r(G^,x). 

8. Finding barely sufficient distributions 

The following result shows how a superset of B{Guj,x) can be constructed using recursion 
starting at Bo{G^,x) = {l{x}}. 

Proposition 8.1. If p G Bi-f-i{Guj,x) then p = Qr-i for some q € Bi{G^j,x) and pebbling 
move r. 

Proof. Suppose that p G Bi+i{Guj,x). Then there is an executable sequence s = 
(si, . . . ,Sj+i) of pebbling moves such that Ps{x) > 1. Then with q = ps^ we clearly have 
p = Qg-T^- Vertex x is reachable from q in the i moves of the sequence (s2, • • • , Sj+i). If x 
is reachable from q in j moves then it is reachable from p in j + 1 moves. So x cannot be 
reached from q in fewer than i moves, which means that q G Bi{Goj,x). □ 

We do not have to use every pebbling move r during the construction of Bi+i^G^j, x) from 
Bi{Guj,x) as shown in the next result that essentially a simple case of the No-Cycle Lemma. 

Proposition 8.2. Let p £ B{G^,x). If ps > o-'nd (v^u) G S then q = p{y^u)-\ 

B{G^,x). 

Proof. Let q{u, v, *) = {q{u) — 1, q{v) — 1, q{*)) and R = S \ {{v^u)}. Then q < q and 
qR{u, V, *) = {qs{u) - 1, qs{v) + 2, qs{*)) 
= iQsiu) - 2, qsiv) + 1, qsi*)) 
= {Psiu),Ps{v),Ps{*)) =Psiu,v,*) 

which means (Jr = ps ^ ^{x}- So q is not barely sufficient. □ 
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Input: Gu), X 
Output: B{Guj,x) 



1 {p,E,W):=m,y,E{G),(l>) I 

2 (5.pushBack((p,£;, VF)) 

3 for (p, E,W) do 
for u € supp(p) do 

for (u, u) £ E and u ^ W do 
Q ■= 

F:=E\{iu,v)} 
X ■.= WU {u} 
for {q, F,X)eQ do 
if g < g then 
|_ break 

ii q = q then 

F := FnF 
X := xnx 
break 

if q > q then 

|_ Q.remove({q,F,X)) 

if did not break then 
|_ (5.pushBack((g,F,X)) 



21 e(G^,rE) :={p| (p,F,W^) gQ} 



(distribution, transfers, forbidden vertices) 
/ / growing queue of distributions 

/ / u has a pebble 
/ / allowed transfer from v to u 
1 1 candidate distribution 
/ / backward transfer no longer allowed 
/ / transfer to u no longer allowed 

/ / candidate too large? 
/ / candidate fails 

/ / candidate already in queue? 
/ / fewer allowed edges for q 
j j not initial in any way 

// g is not barely sufficient? 
/ / remove q from queue 

/ / candidate works 
/ / modification (see Note 19.311 



Figure 8.1. Algorithm to find the set B{Guj,x) of barely sufficient distributions. 



An important interpretation of this result is that every distribution in B{Guj,x) can be 
gotten as pr where p = l^x} and T is a multiset of inverse pebbling moves such that 
{v^u)~^ and {u^v)~^ are not in T together for any u and v. Keeping track of the directions 
of the inverse pebbling moves speeds up the calculation of finding B{G^^x). It also helps 
eliminating moves that cannot be initial moves. 

Given a graph G, let G be the directed graph whose vertex set is V{G) and whose arrow 
set contains two arrows {u,v) and {v^u) for every edge uv € E{G). 

Algorithm 8.3. The algorithm shown in Figure \8l\ finds the set of barely sufficient distri- 
butions. 

The heart of the algorithm is Proposition l8.ll We apply inverse pebbling moves to transfer 
pebbles in barely sufficient distributions in hope of finding new barely sufficient distributions. 
We use triples of the form {p, E, W) where p is a pebble distribution. The role of E is to keep 
track of the direction of the pebble flow so that we can avoid the back and forth transfer 
as explained in Proposition 18.21 The role of W is to avoid pebbling sequences that are 
not regular as explained in Theorem 13.41 (4). Now we give the detailed explanation of the 
algorithm: 
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— lines 1-2: We fill the queue Q of barely sufficient distribution candidates with 
Bq{G^,x). We set E = E{G) since the pebbles can flow in any direction. We 
set = since no vertex is ruled out as the starting vertex of an initial move. 

— line 3: This loop takes an element p oi Q and applies a possible inverse pebbling 
move to create a new distribution q. If p € Bi{x) then g is a candidate for Bi+i{x). 

— line 4: We find a vertex u that has at least one pebble. We plan to remove a pebble 
from this vertex and add two pebbles to an adjacent vertex v. 

— line 5: We only want to apply {v^u)~^ if {u^v)~^ was not used before and if (v->n) 
is an initial move. 

— line 6: We apply the inverse pebbling move {v^u)~^ to create the new barely suffi- 
cient candidate q. 

— line 7: According to Proposition 18.21 we do not want to apply {u-fv)~^ since we 
already used {v^u)~^. 

— line 8: Any move of the form {u^w) is not an initial move since we already have a 
move of the form (v^u). 

— line 9: The loop checks the newly created candidate against the other distributions 
in the queue. 

— lines 10-11: We already put a smaller candidate is the queue so the new candidate 
cannot be barely sufficient. 

— line 12: The new candidate q is already in the queue. It is likely that it was created 
using different inverse pebbling moves. We do not add this candidate to the queue 
twice. Still, we can update the information about this distribution in the queue. 

— line 13: We can reduce the possible inverse pebbling moves using Proposition 18.21 

— line 14: It is possible that a move is initial in one set of pebbling moves but not in 
another set. We only want to declare a move not initial if it is not initial in every 
possible set of pebbling moves that reaches the goal vertex. 

— lines 16-17: If the new candidate is smaller than a distribution q in the queue then 
q cannot be barely sufficient. Therefore we remove it from the queue. 

— lines 18-19: The new candidate is added to the queue. 

Example 8.4. Let x be the goal vertex and vi and V2 be the other vertices of the complete 
graph G = and let uj{xv2) = 5. Figure [8?2l shows how Algorithm 18.31 finds B{G^^x) = 
|p(o)^p(ii)^p(2i)^p(3i) 1 jvjQj^g lYioX p*^^^) is added to the queue and only removed later when 
is found. This late recognition of the fact that p^^^^ is not barely sufficient is the reason 
why the algorithm needs to test p^^^) as a candidate. 

9. Squished distributions 

In this section we prove a version of the Squishing Lemma of [1] using open ears. A pebble 
distribution is squished on a path P if all the pebbles on P are placed on a single vertex of 
P or on two adjacent vertices of P. A pebble distribution can be made squished on a path 
using squishing moves. 

Lemma 9.1. (Squishing) If vertex x is not reachable from a pebble distribution p with size 
n, then there is a pebble distribution of size n that is squished on each unweighted open ear 
and from which x is still not reachable. 

Proof. Let E be an unweighted open ear with consecutive vertices fo, • • • ,Wn- Suppose that 
the pebble distribution p is not squished on E. Let i be the smallest and j be the largest 
index for which p{vi) > and p{vj) > 0. Note that we must have j — i > 2. Define a new 
pebble distribution q by applying a squishing move that moves one pebble from Vi to ffc and 
another pebble from Vj to Vk for some k satisfying i < k < j. 
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1. g = {pW} 

2. Q = {p(0),p(ii)} 

3. Q = {p(0),p("),p(i2)} 

4. Q = {p(0)^p(ll),p(12)^p(21)| 

5. p*^^^-' tested but is larger than p^^^^ 



6. p^^^^ tested, smaller than p^^"^^ 

7. Q = {pW,p(ii),p(2i)} 

8. Q = {p(0),p(ii),p(2i),p(3i)} 

9. p^^2) tested but is larger than p^^^ 

10. p^^^) tested but is larger than p^^^ 



Figure 8.2. The distributions in solid frames belong to B{G^,x). A distri- 
bution p^*-'^ in a dotted frame is never in the queue. A solid arrow from q to 
p is drawn with label (v^u) if q = p(„_>„)-i. A dashed arrow from g to p is 
drawn if g > p and so q ^ B{Guj,x). The table shows how the distribution 
queue Q changes during the execution of the Algorithm 18.31 

Suppose X is reachable from q, that is, there is an acyclic multiset S of pebbling moves 
such that qs > Pick a maximal directed path of T{G,S) with consecutive vertices 

Vk = Wo, wi, ... ,wi all in the set {vi,Vi+i, . . . ,Vj}. Let D be the set of moves corresponding 
to the arrows of this directed path, that is, D = {(wq^wi), . . . , (wi-i^wi)} and let R = 
S \ D. We need to consider three cases depending on whether wi = Vk, wi G {vi.,Vj} or 
wi {vk,Vi,Vj} . It is easy to see that in all three cases we must have pR > Ij^} which is a 
contradiction. 

Applying squishing moves repeatedly on E makes the pebble distribution squished on E. 
This procedure keeps the goal vertex x unreachable. A squishing move on E might remove 
a pebble from another open ear but it cannot add a pebble to it. So if the distribution 
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is squished on an open ear then it remains squished after the application of a squishing 
move on E. So the desired pebble distribution can be reached by applying all the available 
squishing moves on all unweighted ears in any order. □ 

The set 2s{Guj,x) of squished insufficient distributions is the set of those elements of 
I{Guj,x) that are squished on all open ears of G. The set Bs{Guj,x) of squished barely 
sufficient distributions is the set of those elements of B{Goj,x) that are squished on all open 
ears of G. 

Proposition 9.2. We have p £ Is{Guj,x) if and only if q < p does not hold for any 
q G BsiG^,x). 

Proof. The result follows from Proposition 17.31 □ 

Note 9.3. We can find Bs{Gi^,x) by a slight modification of Algorithm 18.31 On line 20 we 
simply remove {p, E, W) from Q if p is not squished. 

Corollary 9.4. Tr{Guj,x) = max{|p| : p G ls{Guj,x)} + 1. 

Proof. The result follows from the Squishing Lemma using tt{Guj,x) = max{|p| : p G 

10. Finding insufficient distributions 

Now we present an algorithm for finding Ts{x). Let p be the pebble distribution defined 
by p{v) = max.{q{v) \ q £ Bs{x)} for all v G V{G). It is clear that if p G Ts{x) then p < m. 
The idea of the algorithm is to decrease the number of pebbles at certain vertices of p until 
it becomes insufficient. 

Algorithm 10.1. The algorithm shown in Figure [Ton\ finds tt{Guj,x) using Bs{Guj,x). 

The algorithm uses Proposition 19.21 and Corollary 19.41 The input Bs{Guj,x) is the output 
of the modified Algorithm 18.31 as explained in Note 19.31 It contains all the squished barely 
sufficient distributions. Now we give the detailed explanation of the algorithm: 

— lines 1-3: For each vertex v there is a barely sufficient distribution that has pebbles 
only on v. So the distribution p is an upper bound for T{Gu},x). Every insufficient 
distribution can be constructed form p by decreasing the number of pebbles on some 
vertices. 

— lines 4-12: We only need to find those insufficient distributions that are squished. So 
instead of using p we can use a new distribution p by removing all the pebbles from p 
at a few vertices until the distribution becomes squished. Every squished insufficient 
distribution can be constructed from one such p by decreasing the number of pebbles 
on some vertices. So we collect all these p's in the queue Q on line 7. The queues P 
and Q are kept ordered so lines 7 and 12 can use binary searches to avoid duplication 
in the queues. 

— line 13: The variable M is initialized, it will contain the maximum size max{|p| : 
p G Is{Guj,x)} of an insufficient distribution. 

— lines 14-27: We are finding M using Corollaries 19.21 and 19.41 Queue Q contains 
pairs of the form {p,i). Such a p is always smaller than the first i — 1 elements of 
Bs{Guj,x). In each iteration we replace {p,i) G Q by possibly several new elements 
in the queue of the form {q, i + 1). Such a g is now smaller than the first i elements 
of Bs{Guj,x). We can find such a g by decreasing the number of pebbles at a vertex 
where the i-th element of Bs{Guj,x) has pebbles. The loop in line 23 finds these 
vertices. Line 27 uses binary search. 

— line 28: The pebbling number is calculated according to Corollary 19.41 
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Input: C := Bs{Gi_j 
Output: tt{G^,x) 



1 for V G V{G) do 

2 |_ p{v) := Ta.ayi{q{v) \ q ^ C} 

3 P.pushBack(p) / , 

4 while P not empty do 
P.popBack(p) 
if p squished then 

Q. insert ((p, 1)) 
continue 

for V e V{G) do 
q:=p 
q{v) := 
P. insert iq) 

13 M := 

14 while Q not empty do 
(5.popBack((p, i)) 
if IpI < M then 

|_ continue 

while i < \C\ and C[i] ^ p do 
\_i:=i + l 

if i > \C\ then 
M := IpI 
continue 

for V G V{G) do 
g :=p 

:=C[i](t;)-l 
if q{v) > then 
|_ Q.iiisert({q,i + 1)) 

28 tt{G^,x) := M + 1 



// p is an upper bound for Is{Gi^,x) 

ordered queue of not yet squished candidates 
/ / more candidates to try 
/ / work with candidate p 

1 1 ordered queue of squished candidates 
/ / nothing more to do with p 

1 1 try to improve p 
1 1 modify p 
1 1 this might make it squished 
/ / add improved candidate to the queue 

// size of best insufficient distribution so far 
/ / more candidates to try 
1 1 p works for . . . ,C[z — 1] 
// too few pebbles? 
/ / candidate has no hope to be better 

/ / find first i such that C[i] <p 
II not found yet 

//no such i, candidate works 
// p is the best insufficient distribution so far 
/ / nothing more to do with p 

1 1 make candidate work for C[z] 
/ / modify p 
1 1 q works for C [1] , . . . , C [i] 
/ / nonnegative number of pebbles on t>? 
/ / add improved candidate to the queue 



Figure 10.1. Algorithm to find the distribution with the most 
pebbles that is insufficient for the goal vertex. 



It is important to keep our queues sorted and to use binary search at the insert operations. 
Without this the algorithm becomes too slow to be practical. 



11. Test results 

We tested our algorithms by calculating the pebbling number of every connected graph 
with fewer than 10 vertices. We used Nauty [7j to generate these graphs and their auto- 
morphism groups. We simplified each graph as follows. For each goal vertex we replaced 
every closed ears and cut ears with a weighted edge as described in Subsections 16.11 and 
16. 2[ Then we recursively used available leaves to simplify the graph as much as possible as 
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\V{G)\ = 1 

1 1 

\V{G)\=2 

2 1 

\V{G)\=3 

3 1 
4 1 

|F(G)|=4 

4 3 

5 2 



|y(G)| = 5 

i 10 
i 5 




\V{G)\ = 7 
7 322 







9 


1658 


8 


113 


10 


1870 


9 


125 


11 


1425 


10 


129 


12 


478 


11 


68 


13 


26 


12 


4 


14 


1 


16 


23 


16 


190 


17 


35 


17 


341 


18 


22 


18 


333 


19 


2 


19 


148 






20 


15 


32 


4 






33 


5 


32 


36 






33 


52 


64 


1 


34 


34 


\V{G)\=8 


35 


3 


8 


4494 


64 


6 



65 



128 



\V{G)\ 



9 126646 

10 43935 

11 41222 

12 22756 

13 4975 



14 
15 

16 
17 
18 
19 
20 
21 



208 
6 

2505 
5293 
5992 
4070 
1310 
137 



22 
23 

32 
33 
34 
35 
36 
39 

64 
65 
66 
67 

128 
129 

256 



Table 1. The frequency of pebbling numbers for graph with less than 10 
vertices. The data is grouped by the number of vertices in the graph. 

Each data row contains a possible pebbling number followed by the frequency 
of this pebbling number. 



described in Section [5l Then we run Algorithms 18.31 and 110.11 to find the pebbling number 
of the simplified graph. 

The automorphism group helped us reducing the number of goal vertices to representatives 
of orbits. It is well known that the hardest to reach goal vertex in a tree is a leaf, so in trees 
we only picked leaf vertices for the goal vertex. 

The algorithms were coded in C++ using the Standard Template Library. The code was 
compiled with the gnu compiler. It took about a day on a 3 GHz Unix machine to finish the 
calculations. The calculation for graphs with fewer than 9 vertices took less than 10 minutes. 
Table [U shows the frequency of the pebbling numbers. The result confirms the existence of 
gaps in the spectrum of pebbling numbers described in [2]. We checked our results on many 
graphs with known pebbling numbers such as paths, complete graphs, cycles, some trees, 
Lemke graph, Petersen graph. The pebbling numbers are available on the internet. 

12. Further questions 

(1) We do not have any example where t i— > 7rt{G,x) is not linear for t > 3. Is it true 
that t I— > 7rt{G,x) is always linear for t > to for some to? What properties of the 
graph can be used to find to? 

(2) Is it possible that the goal vertex x that maximizes n{G, x) is an interior vertex of 
a cut ear of G? The answer seems to be no and it may depend on the first question. 
We could use this result to speed up our calculation of 7r(G) since we would have to 
test fewer goal vertices. 

(3) Is it possible to take advantage of open ears in a better way? Can we simplify a 
graph with an open ear so that it has fewer vertices? This simplification might help 
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tremendously if it is compatible with ear decomposition. It might be possible to 
reduce a graph completely if we know tt^ for all graphs with fewer vertices. 

(4) Perhaps adding extra weighted edges could speed up Algorithm 18.31 in certain cases. 
For example in an open ear, we could connect the end vertices to the interior vertices 
with appropriate weights depending on the distance. 

(5) What general results are there about the pebbling number of weighted graphs? In 
particular, does Graham's conjecture hold for weighted graphs? 

(6) What is the pebbling number of simple weighted graphs like the weighted complete 
graph or a weighted cycle? Even a weighted triangle seems to be fairly complicated 
with a lot of cases to consider. 
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